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O ■ ABSTRACT 

■ In the present lore of cosmology, matter and space-time emerged from a singularity and 
I ' 

Q . evolved through four different regimes: inflation, radiation, dark matter and dark energy 
H ■ 

c/3 ' dominated eras. In the radiation and dark matter dominated stages, the expansion of the 

: 

Universe decelerates while the inflation and dark energy eras are accelerating regimes. So 
far there is no clear cut connection between these accelerating periods. More intriguing, the 
substance driving the present accelerating stage remains a mystery, and the best available 
candidate (A-vacuum) is plagued with the coincidence and cosmological constant problems. 
^ \ In this paper we overcome such problems through an alternative cosmic scenario based 

on gravitationally-induced particle production. The model proposed here is non-singular 
with the space-time emerging from a pure initial de Sitter stage thereby providing a nat- 
^ ' ural solution to the horizon problem. Subsequently, due to an instability provoked by the 
production of massless particles, the Universe evolves smoothly to the standard radiation 
dominated era thereby ending the production of radiation as required by the conformal in- 
variance (Parker's theorem). Next, the radiation becomes subdominant with the Universe 
entering in the cold dark matter dominated era. Finally, the negative pressure associated 
with the creation of cold dark matter particles accelerates the expansion and drives the 
Universe to a final de Sitter stage. The late time cosmic expansion history is exactly like in 
the standard ACDM model, however, there is no dark energy. The model evolves between 
two limiting (early and late time) de Sitter regimes. Our scenario is fully determined by two 
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extreme energy densities, or equivalently, the associated de Sitter Hubble scales connected 
by p:/pf = {H,/Hff ~ 
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The microscopic description for gravitationally-induced particle production in an expand- 
ing Universe began with Schrodinger's [l| seminal paper, who referred to it as an alarming 
phenomenon. In the late 1960s, this issue was rediscussed by Parker and others 2, ^ based 
on the Bogoliubov mode-mixing technique in the context of quantum field theory in curved 
space-time. Physically, one may think that the (classical) time varying gravitational field 
works like a 'pump' supplying energy to the quantum fields. 

In order to understand the basic approach, let us now consider a real minimally cou- 
pled massive scalar field evolving in a flat expanding Friedman- Robert son- Walker (FRW) 
geometry. The field is described by the following action^ 



S = - I yj—gd'^x 



In terms of the conformal time r] {dt = a{ri)dri), the metric tensor g^^ is conformally equiv- 
alent to the Minkowski metric rj^^,, so that the line element is ds"^ = a^{ri)ri^jjdx^dx'^ , where 
a{ri) is the cosmological scale factor. Writing the field (j){r],x) = a{ri)~^Xj obtains from 
the above action 

x"-V2x+(mV-^)x = 0, (2) 

where the prime denotes derivatives with respect to rj. Notice that the field x obeys the 
same equation of motion as a massive scalar field in Minkowski space-time, but now with a 
time dependent effective mass, 

a" 

mljf{ri) = m^a^ - — . (3) 

This time varying mass accounts for the interaction between the scalar and the gravitational 
fields. The energy of the field x is not conserved (its action is explicitly time-dependent), 
and, more important, its quantization leads to particle creation at the expense of the classical 
gravitational background 2-^. 

On the other hand, in the framework of general relativity theory (TRG), the scale factor 
of a FRW type Universe dominated by radiation satisfies the following equation: 



ad + d = (4) 



We adopt units such that h — ks — c — 1. 
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or, in the conformal time, a" = 0. Therefore, for massless fields (m = 0), there is no particle 
production since Eq. (|2]) reduces to the same of a massless field in Minkowski spacetime, 
and, as such, its quantization becomes trivial. This is the basis of Parker theorem concerning 
the absence of massless particle production in the early stages of the Universe. Note that 
Parker's result does not forbid the production of massless particles in a very early de Sitter 
stage (a" 7^ 0). Potentially, we also see that massive particles can always be produced by 
a time varying gravitational field. As we shall see, such features are incorporated in the 
scenario proposed here. 

In principle, for applications in cosmology, the above semiclassical results has three basic 
difficulties, namely: 

(i) The scalar field was treated as a test field, and, therefore, the FRW background is not 
modified by the newly produced particles. 

(ii) The particle production is an irreversible process, and, as such, it should be con- 
strained by the second law of thermodynamics. 

(iii) There is no a clear prescription of how an irreversible mechanism of quantum origin 
can be incorporated in the Einstein Field Equations (EFE). 

Later on, a possible macroscopic solution for these problems was put forward by Prigogine 
a„d coworkers Q using no,r-equ,libr,um thernrody,ranr,cs for open systems, and by Calvao, 
Lima & Waga [6| through a covariant relativistic treatment for imperfect fiuids (see also [7[). 
The leitmotiv of the approach is that particle production, at the expense of the gravitational 
field, is an irreversible process constrained by the usual requirements of non-equilibrium 
thermodynamics. This irreversible process is described by a negative pressure term in the 
stress tensor whose form is constrained by the second law of thermodynamics^. 

In comparison to the standard equilibrium equations, the irreversible creation process is 
described by two new ingredients: a balance equation for the particle number density and 
a negative pressure term in the stress tensor. Such quantities are related to each other in a 
very definite way by the second law of thermodynamics. Since the middle of the nineties, 
several interesting features of cosmologies with creation of cold dark matter and radiation 



have been investigated by many authors 



Hul. 



^ The semiclassical approach is unable to provide the entropy burst accompanying the particle production 
since it is adiabatic and reversible. 
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In this context, we are proposing here a new cosmological scenario where the accelerating 
stages of the cosmic evolution are powered uniquely by the creation of massless and massive 
cold dark matter particles. In this model, the Universe starts from a de Sitter dominated 
phase (a oc e^'^) powered by the production of massless particles. Subsequently, it deflates 
and evolves to the standard radiation phase (a oc t^^^) thereby ending the creation of massless 
particles. Due to expansion, the radiation becomes subdominant with the Universe entering 
in the cold dark matter (CDM) dominated era. Finally, the negative pressure associated 
with the creation of cold dark matter particles accelerates the expansion and drives the 
Universe to a final de Sitter stage. The horizon problem is naturally solved in the initial 
de Sitter phase. In addition, the transition from Einstein-de Sitter (a oc 

^2/3) to a de Sitter 

final stage (a oc e^^^) guarantee the consistence of the model with the supernovae type la 
data and complementary observations. A transition redshift of the order of a few (exactly 
the same value predicted by ACDM) is also obtained. 

For simplicity, let us consider the EFE for a fiat geometry: 

• 2 •■ • 2 

SnGp = 3^, SnGip + p,) = -2- - ^, (5) 

where an overdot means time derivative, p and p are the dominant energy density and 
pressure of the cosmic fiuid, respectively, and Pc is a dynamic pressure which depends on 
the particle production rate. Special attention has been paid to the simpler process termed 
"adiabatic" particle production. It means that particles and entropy are produced in the 
space-time, but the specific entropy (per particle), a = S/N, remains constant [6']. In this 
case, the creation pressure reads [6-9] 

Pc = (6) 

where F with dimensions of {time)~^ is the particle production rate and H = a/a is the 
Hubble parameter. 

How the evolution of a{t) is affected hyV? By assuming a dominant cosmic fiuid satisfying 
the equation of state (EoS), p = up, where w is a constant, the EFE imply that 

H+l{l + oo)H^(l-^)=0. (7) 



2' ' \ 3HJ 

The de Sitter solution {H = 0, T = 3H = constant) is now possible regardless of the EoS 
defining the cosmic fiuid. Since the Universe is evolving, such a solution is unstable, and, 
as long as F << 3H, conventional solutions without particle production are recovered. 



=^ The main effect of T is to provoke a dynamic instability in the space-time, thereby 
allowing a transition from de Sitter to a conventional solution, and vice versa. 
A. Prom an early de Sitter stage to the standard radiation phase 

Let us first discuss the transition from an initial de Sitter stage to the standard radiation 
phase. The main theoretical constraints are: 

• The model must not only solve the horizon problem but also provide a quasiclassical 
boundary condition to quantum cosmology ( a hint on how to solve the initial singularity 
problem). 

• Massless particles cannot be quantum-mechanically produced in the conventional radi- 
ation phase (Parker's Theorem). 

To begin with, let us assume a radiation dominated Universe {u = 1/3, F = F^). The 
dynamics is determined by the ratio Tr/SH (see ([7])). The most natural choice would be a 
ratio which favors no epoch in the evolution of the Universe (F,,/3if = constant). However, 
the particle production must be strongly suppressed, Tr/3H << 1, when the Universe enters 
the radiation phase. This means that the expansion The simplest formula satisfying such a 
criterion is linear, namely: Tr/SH = H/Hj, where Hj is the initial de Sitter expansion rate 
{H < Hi). Inserting this into ([7]) it becomes: 

H + 2H^ (^1 - -^^ = 0. (8) 
The solution of the above equation can be written as 

where D > is an integration constant. Note that if = if/ is a special solution of Eq.(|8]) 
describing the exponentially expanding de Sitter space-time. This solution is unstable with 
respect to the critical value D = 0. For D > 0, the universe starts without a singularity and 
evolves continuously towards a radiation stage, a ~ t^/^, when Da^ >> 1. By integrating 
we obtain the scale factor: 

Hit = In - + ^{a/a,)\ (10) 
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where = Dal is an integration constant and defines the transition from the de Sitter 
stage to the beginning of the standard radiation epoch^. At early times (a ^ a*), when the 
logarithmic term dominates, one finds a ~ a^,e^^*, while at late times, a ^ a^,, H Hj, 
flTU]) reduces to a ^ a* (^^t) ^ , and the standard radiation phase is reached. 



It should be noticed that the time scale Hj ^ provides the greatest value of the energy 
density, pj = characterizing the initial de Sitter stage which is supported by the maximal 
radiation production rate, = 3Hj. From ([5]) and (Q we obtain the radiation energy 
density: 



Pr = Pi 



1 + ( - 

a* 



1 -2 



(11) 



As expected, we see again that the conventional radiation phase, pr ^ a ^, is attained when 
a a*. 

How the cosmic temperature evolves? For "adiabatic" particle production the energy 
density scales as pr ~ T'^ js, 7|, and the above equation implies that 



Tr = Tj 



-1/2 



(12) 



where T/ is the temperature of the initial de Sitter phase which must be uniquely determined 
by the scale Hj. We see that the expansion proceeds isothermally during the de Sitter phase 
(a << a*). A basic consequences is: 

=^ The supercooling and subsequent reheating taking place in several inflationary variants 
are avoided In other words, there is no the so-called 'graceful exit' problem. 

After de Sitter stage, the temperature decreases continuously in the course of the expan- 
sion. For a » a^: {H << Hj), we obtain T ~ a~^. Accordingly, the comoving number of 
photons becomes constant since n oc a^^, as expected for the standard radiation stage. ^ 

What about the initial temperature Tj ? Since the model starts as a de Sitter space- 



time, the most natural choice is to define Tj as the Gibbons-Hawking temperature 



isj of 



^ This kind of evolution was first discussed by G. L. Murpliy [12| by studying possible effects of the second 
viscosity in the very early Universe. Later on, it was also investigated in a more general framework 



involving cosmic strings by J. D. Barrow 13| who coined the exp ression "Deflationary Universes" . It has 
also been discussed in connection with decaying A(i)-models 14 1. 

3 
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Since rir (x T^, the average photon concentration reads rir = nj 



1 + A- 
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its event horizon, Tj = Hi/2'k. Naively, one may expect T/ of the same order or smaller 
than the Planck temperature because of the classical description. From EFE we have pi = 
Smpi'^Hj'^ /8n (where mp; ~ 1.22 x lO^^GeV), and since the energy density is p/ = N^(T)Tj, 
one finds Tj Hj lO^^GeV (where A^*(T) = 7r^(7*(T)/30 depends on the number of 
effectively massless particles). 

Naturally, due to the initial de Sitter phase, the model is free of particle horizons. A 
light pulse beginning at t = — oo will have traveled by the cosmic time t a physical distance, 
dnit) = Ci{t) j-oo which diverges thereby implying the absence of particle horizons: 

^ The local interactions may homogenize the whole Universe. 

Since photons are not produced in the radiation phase, the Big-Bang Nucleosynthesis 
(BBN) may work in the conventional way 
dark matter (Einstein-de Sitter, a{t) oc t^/^) dominated phase. 

B. Prom Einstein-de Sitter to a late time de Sitter stage 

Due to conservation of baryon number the remaining question is the production rate of 
cold dark matter particles and the overall late time evolution. In other words, what is the 
form of TdrrJ For simplicity, we consider here only the dominant CDM component. 

In principle, should be determined from quantum field theory in curved spacetimes. 
le absence of a rigorous treatment, we consider (phenomenologically) the following fact 



n t 



19 



All available observations are in accordance with the ACDM evolution both at the 
background and perturbative levels. 



Now, we recall that a flat ACDM model evolves like: 



1 - 



H 



0, (13) 



where Hf^ = A/3 sets the Hubble scale of the final de Sitter stage {H > Hj). Such behavior 
should be compared to that predicted for a dust filled model [u = 0, T = T^m) with particle 
production (see Eq. ([7])): 

H + (l - ^) = 0. (14) 

By comparing ( fT3l) and ( fT4l) . we see that the same background evolution requires that 
T dm I'm = [Hf/Hf'. The limiting value of the creation rate, Tdm = '^Hf, leads to a late 
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FIG. 1: (a) The likelihood for 17^ based on 307 Supernova data (Union) [l^. (b) Evolution of 
the scale factor predicted by the matter creation model (solid line) and the traditional ACDM 
cosmology (open points). In this plot we have adopted the best fit, Cl\ = 0.72, from Supernova 
data. 

time de Sitter phase {H = 0, H = Hf) thereby showing that the de Sitter solution now 
becomes an attractor at late times. With this proviso, the solution of (IT^ reads: 



2 



:i-nj,){i + zy + nA\, (15) 

where Cl\ = [Hf/HoY is smaller than unity and 1 + z = a~^. Such solution mimics the 
Hubble function H{z) of the traditional fiat A-cosmology, with VL\ playing the dynamical 
role of VLf^ (dark energy appearing in the concordance model). The dark matter parameter 
{^dm = 1) is also replaced by an effective parameter, {Qdm)eff = 1 — ^a, which quantifies 
the amount of matter that is clustering. This explains why this model is in agreement with 
the dynamical determinations related to the amount of the cold dark matter at the cluster 
scale, and, simultaneously, may also be compatible with the position of the first acoustic 
peak in the pattern of CMB anisotropics which requires Qtotai = 1- 
By integrating (llSp we obtain: 

1/3 



=^ The late time dynamics is determined by a single parameter (Qa) and is identical to 
that predicted by the flat ACDM model. 
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In Figure 1, we show the hkehhood of Qa based on the Union supernova sample. Note 
also that by replacing the value of Tdm into the definition of the creation pressure (see Eq. [6]) 
one obtains that it is negative and constant {pc = —SHj/SnG = —SQaHq/SttG). Therefore, 
the late time evolution of our c omp lete cosmological scenario coincides exactly with the one 



recently discussed in Refs. 



22 



23| following a slightly different approach. 



Concluding, a new cosmology based on the production of massless particles (in the early 
de Sitter phase) and CDM particles (in the transition to a late time de Sitter stage) has 
been discussed. The same mechanism avoids the initial singularity, particle horizon and the 
late time coincidence problem of the ACDM model has been eliminated (A = 0). 

In this scenario, the standard cosmic phases - a radiation era followed by an Einstein-de 
Sitter evolution driven by nonrelativistic matter until redshifts of the order of a few - are not 
modified. However, the model has two extreme accelerating phases (very early and late time 
de Sitter phases) powered by the same mechanism (particle creation). Therefore, it sheds 
some light on a possible connection among the different accelerating stages of the universe. 
In particular, since Hj = Q^Hq, where ~ 0.7 and Hq — 1.5 x 10~'^'^GeV, it sets the ratio 
of the primeval and late time de Sitter scales to be pi/pf = {Hj/Hf)"^ ~ 10^^^. Such a result 
in the present context has no correlation with the so-called cosmological constant problem 



1251. 

As it appears, the cosmic history discussed here is semi-classically complete. However, 
there is no guarantee that the initial de Sitter configuration is not only the boundary con- 
dition of a true quantum gravitational effect. In other words, the very early de Sitter phase 
may be the result of a quantum fluctuation which is further semi-classically supported by 
the creation of massless particles (in this connection see 2J] and Refs. there in). 

Naturally, the existence of an early isothermal de Sitter phase suggests that thermal 
fluctuations (within the de Sitter event horizon) may be the causal origin of the primeval 
seeds that will form the galaxies. Such a possibility and its consequences for the structure 
formation problem deserves a closer investigation and is clearly out of the scope of the 
present paper. 

At present, we also know that a more complete version of the late time evolution must be 
filled with CDM (~ 96%) and baryons (~ 4%), and, unlike A-cosmology, the baryon to dark 



matter ratio is a redshift function 



22 



231]. In particular, this means that studies involving 



the gas mass fraction may provide a crucial test of our scenario, potentially, modifying our 
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present view of the dark sector. Some investigations along the above discussed hnes are in 
progress and will be published elsewhere. 
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